Bulk Viscous Matter-dominated Universes: Asymptotic Properties 
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By means of a combined study of the type la supernovae test,together with a study of the asymp- 
totic properties in the equivalent phase space - through the use of the dynamical systems tools - we 
demonstrate that the bulk viscous matter-dominated scenario is not a good model to explain the 
accepted cosmological paradigm, at least, under the parametrization of bulk viscosity considered in 
this paper. The main objection against such scenarios is the absence of conventional radiation and 
matter-dominated critical points in the phase space of the model. This entails that radiation and 
matter dominance are not generic solutions of the cosmological equations, so that these stages can 
be implemented only by means of very particular solutions. Such a behavior is in marked contra- 
diction with the accepted cosmological paradigm which requires of an earlier stage dominated by 
relativistic species, followed by a period of conventional non-relativistic matter domination, during 
which the cosmic structure we see was formed. 

PACS numbers: 04.20.Ha, 98.80.-q, 98.80.Es, 98.80.Jk 



I. INTRODUCTION 



Within the context of early inflation, it has been known 
since long time ago that an imperfect fluid with bulk vis- 
cosity in cosmology can produce an accelerated expan- 
sion without the need of a cosmological constant or some 
other inflationary scalar field [l| (although some authors 
do not agree with this conclusion 0]). Further extrapola- 
tion of this idea - used to induce an accelerated expand- 
ing Universe without the need of unknown components 
- leads to the possibility that one alternative candidate 
to explain the present acceleration can be bulk viscous 
pressure No matter how attractive it seems, this idea 
faces some problems, among them, the need to have a 
satisfactory mechanism for the origin and composition of 
the bulk viscosity (see, for instance [1, [1]). 

From a thermodynamical point of view the bulk vis- 
cosity in a physical system is due to its deviations from 
the local thermodynamic equilibrium |5[ . In a cosmologi- 
cal setting, the bulk viscosity may arise when the cosmic 
fluid expands (or contracts) too fast so that the system 
does not have enough time to restore its local thermody- 
namic equilibrium and then it arises an effective pressure 
restoring the system to its thermal equilibrium. When 
the fluid reaches again the thermal equilibrium then the 
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bulk viscous pressure vanishes |5, 6]. Therefore, in an 
accelerated expanding Universe, it may be natural to as- 
sume the possibility that the expansion process is actu- 
ally a collection of states out of thermal equilibrium in 
a small fraction of time giving rise to the existence of a 
bulk viscosity Q. 

Usually the way to test the theoretical (and observa- 
tional) viability of a given cosmological model is through 
using known solutions of the cosmological field equations, 
or by seeking for new particular solutions that are phys- 
ically plausible. However, in general, the cosmological 
field equations are very difficult to solve exactly. Even 
if a given analytic solution can be found, it will not be 
unique, but just one in a large set of them. This is not 
to talk about stability of given solution(s). An alterna- 
tive way around is to invoke the dynamical systems tools 
to extract very useful information about the asymptotic 
properties of the model instead. In this regard, knowl- 
edge of the critical (also equilibrium or singular) points in 
the phase space - corresponding to a given cosmological 
model - is a very important information since, indepen- 
dent on the initial conditions chosen, the orbits of the 
corresponding autonomous system of ordinary differen- 
tial equations (ODE) will always evolve for some time in 
the neighborhood of these points. Besides, if the point 
were an (global) attractor, independent of the initial con- 
ditions, the orbits will always be attracted towards it (ei- 
ther into the past or into the future). Going back to the 
original cosmological model, the existence of the equilib- 
rium points can be correlated with generic cosmological 
solutions that might really decide the fate and/or the 
origin of the cosmic evolution. 

A phase space for a model which is consistent with the 
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presently accepted cosmological paradigm should contain 
critical points associated with: i) a radiation-dominated 
(relativistic) stage, followed by ii) a matter-dominated 
(non-relativistic) phase, which is important to allow for 
the formation of the amount of cosmic structure we see, 
and iii) an accelerated expanding stage which possibly 
might last for ever. Perhaps, there should be also a crit- 
ical point in the phase space associated with an early 
inflationary period, however, this would require of fur- 
ther refinements of a given cosmological model which is 
primarily intended to explain the period lasting between 
decoupling of radiation and baryons and up to the present 
accelerating phase. 

In the present paper we will make a combined use of the 
dynamical systems tools and of the type la supernovae 
test, to extract as much as possible useful information 
about the asymptotic properties of a bulk viscous matter- 
dominated Universe, in order to be able to judge about 
its theoretical viability to accommodate the accepted cos- 
mological paradigm. In the specific model we shall be 
investigating the cosmological dynamics is fuelled by a 
conventional matter component jointly with a fluid with 
bulk viscosity of the form: C — Co + (iH + (20/0, where 
Co ; Ci > an d C2 ■ ar e constants to be determined by the ob- 
servations and H = a/a is the Hubble parameter. The 
term £0 takes into account the simplest parametrization 
for the bulk viscosity: a constant. The term (iH charac- 
terizes the possibility of a bulk viscosity proportional to 
the expansion ratio of the Universe, while the third term 
(,2 'd/a takes into consideration the influence acceleration 
of the expansion might has on the bulk viscosity. 

The results of the present study will convincingly show 
that the model of bulk viscosity considered here (see 
also 0,3) is in marked contradiction with the presently 
accepted cosmological paradigm and, hence, should be 
ruled out. In fact, it will be shown that there are not 
any critical points in the phase space of this cosmologi- 
cal model which could be associated with either radiation 
or matter domination. This, in turn, entails that radia- 
tion and matter dominance are not generic solutions of 
the model. These can be, at most, very particular so- 
lutions that can be achieved under specific initial condi- 
tions. This is in marked contradiction with the accepted 
cosmological paradigm in that the latter requires of an 
earlier stage dominated by relativistic species (specifi- 
cally a radiation-dominated era), followed by a period of 
conventional matter domination during which the cosmic 
structure we see was formed. 

The paper has been organized as it follows. The rel- 
evant details of the model are exposed in section [Til fol- 
lowed by its observational testing using the supernovae 
SNe la data in section Ull Al In the first part of the dy- 
namical systems study fsection llV[) . for sake of simplicity, 
we will resort to a single component to embrace radia- 
tion, baryons, etc. Then, in sectionfVl we consider a more 
physically involved model where the cosmic dynamics is 
fuelled by radiation, dark matter, and pressureless bulk 
viscous matter. In order to draw the phase portraits we 



use the best estimated values of the free parameters of 
section Ull Al A thorough discussion of the properties of 
the model in the equivalent phase space is given in sec- 
tion [VT] It will be evident that, independent of the values 
of the free parameters, there are not critical points in the 
phase space that could be correlated with either radia- 
tion or matter dominance stages. In the final section fVIII 
brief conclusions are given. 

II. COSMOLOGY OF BULK VISCOUS 
MATTER-DOMINATED UNIVERSES. 

In this section we analyze a cosmological model com- 
posed by a baryon matter and a bulk viscous components 
with dust behavior (w = 0) as dark matter, in a spatially 
flat universe. The energy-momentum tensor for the pres- 
sureless baryon component is given as a perfect fluid as 
usual 

T (m ) — n 11 11 

J - fj,u — fm 

where p m is the energy density of the matter compo- 
nent and 11^ is the four-velocity vector. The energy- 
momentum tensor of the bulk viscous component is that 
of an imperfect fluid with a first-order deviation from the 
thermodynamic equilibrium. It can be expressed as Q: 

where 

p; = p v - cvuu v , (i) 

p v and P v are the energy density and pressure of the vis- 
cous fluid respectively. The term g^ v is the metric tensor, 
and the subscript V stands for "viscous" component. 
The term £ is a bulk viscous coefficient that arises in 
a fluid when it is out of the local thermodynamic equi- 
librium and which induces a viscous pressure — C^J u u v 
[5|. The term P* is an effective pressure composed by 
the pressure P v of fluid plus the bulk viscous pressure. 
It was initially proposed in reference [10| for relativistic 
dissipative processes in thermodynamic systems out of 
local equilibrium, and later on the authors of Ref. [ll| 
developed an equivalent formulation. 

Here we shall explore the cosmological implications of 
such a model by assuming that the Universe is filled 
with bulk viscous fluid only. For Friedmann-Robertson- 
Walker (FRW) spacetimes with flat spatial sections we 
have 

ds 2 = -dt 2 + a 2 {t){dr 2 + r 2 dVL 2 ), 

where a(t) is the scale factor. The conservation equa- 
tions for the pressureless baryon matter and the viscous 
component can be written respectively as 

Pm + 3Hp m = 0, 

p v + 3H(p v +P:)=0, (2) 
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where H = a I a is the Hubble parameter and, as usual, 
the over-dot stands for time derivative. The solution 
of the first equation in ([2]) for the matter component is 
Pm = PmO a ~ 3 i where p m Q is the present value of the mat- 
ter energy density. On the other hand, the bulk viscous 
pressure — ('V„it 1 ' can be written as —3(,H. So, assuming 
the general equation of state for the viscous component 
as P v — wp v , the second conservation equation in © for 
the viscous fluid becomes 



Co + CiE - (9C 2 /#) 


Sl m oa 3 + (1 + 3w)£l v 


.(. 


-scV) 



(9) 



with C2 7^ 1/9. Hence, the conservation equation for the 
viscous component ([3]) in terms of dimensionless quanti- 
ties is an ordinary differential equation 



p v + 3£T(1 + w) Pv - %H 2 = 0. (3) 
The Einstein's equations for the model arc 



H 2 = 


8irG , 
3 {P 


a 


AttG 




a 


3 


P 



(4) 



We assume the parametrization for the bulk viscosity £ 
of the viscous component p v in the form of the following 
expansion: 



C = Co + G# + C2 



(5) 



where Co? Ci and C2 are constants to be determined by 
the observations. The term (a/a) in ([5]) can be written 
as (a/a) = (d/a)/H. So, using second equation in (@}, it 
can be written as 



4ttG r 



p m + (1 + 3w) Pv - 9(H 



(6) 



Next we substitute eq. © into ([5]) , and after re-arranging 
terms we obtain 



Co + CiH - C2 \p m + (1 + 3w)p v 



1 - 9 



'AttG 



)c 2 



(7) 



Defining the dimensionless bulk viscous coefficients as 



Ci = (24ttG)Ci, £2 = (^)Ca, 

and the dimensionless Friedmann equation - first equa- 
tion in flU - as 



a^+3(l + w)Q v -9(E = Q, (10) 
da 

where E and C are given by the equations © and © 
respectively. Using the relation 1 + z = a -1 between the 
redshift and the scale factor, we express the ODE (fTU|) as 



(l + z)^^--3(l + w)h v (z) + 



9C 



n v ( z ) + n m0 (i + zf 



1/2 



= 0. (11) 



We solve numerically this ODE assuming a dust behav- 
ior for the viscous matter (i.e., w — 0) and the initial 
conditions CI v q = Cl v (z = 0) = 0.96, S! m o = 0.04. 

The dimensionless Hubble parameter E = H/Hq be- 
comes 

E(Co, Ci, C2) = y/fl v (z) + n m0 {l + z)3, (12) 

where Q v (z) is given by the numerical solution of the 
ODE pT)|. with C 2 # 1/9. The dimensionless density 
parameters of the viscous dark matter VL v (z) and baryon 
^m(z) components, can be expressed as 



Q v (z) 
Q m (z) 



n v (z) + n m0 (i + z) 3 ' 

n m0 (i + z) 3 
n v ( z ) + n m0 (i + z) 3 '' 



(13a) 
(13b) 



where (l v (z) is given by the ODE (fTTj) . The evolution of 
ft v (z) and f2 m (z) is shown in figure [SJ For the decelera- 
tion parameter q(a) = — (a/a)H~ 2 1 one obtains 



n m0 (i + z) 3 + a v ( z ) -%e( z ) 
2(O m0 (l + zf + a v (z)) : 



(14) 



where C and E(z) are given by Eqs.© and (JT2]) respec- 
tively. The evolution of q(z) is shown in figure [6l 



as-. 



(8) 



where fi m0 = p m0 /p° rit and 0„ = p v /p" r iv tne equation 
(|7|) can be rewritten in dimensionless form: 



III. COSMOLOGICAL PROBES 

We test the viability of the model and constrain its 
free parameters (Co, Cb C2) using the type la Supernovae 
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FIG. 1: Confidence intervals for the dimensionless viscous co- 
efficients (Co, CO of a bulk viscosity parametrization of the 
form C = Co + (^iH (we have set C2 = 0). We find as best 
estimates (Co = 3.12±°;f|, Ci = -0.96 ± 0.3), see table U The 
solid red, dashed and filled contour plots correspond to the 
use of SNe, H(z) and the joint SNe + H(z) data sets respec- 
tively. The gray shaded area at the bottom left indicates the 
forbidden region where the total bulk viscosity ("total given by 
cq.© is negative at the present-day (z = 0). The confidence 
intervals shown correspond to 68.3%, 95.4% and 99.73% of 
probability. The Hubble constant Ho was marginalized as- 
suming a constant prior probability distribution for Hq. It 
was assumed the values Cl v (z = 0) = 0.96, fi m o = 0.04 and 
w = 0. 



(SNe la) observations and the Hubble parameter H(z) 
measured at different redshifts. We compute the best es- 
timated values for pairs of ((0,(1,(2), the goodness- of- fit 
of the model to the data and the confidence intervals by a 
X 2 function minimization, to constrain their possible val- 
ues with levels of statistical confidence which are shown 
in figures [T] [3] 



A. Type la Supernovae 

We use the "Union2.1" SNe la data set (2012) from 
"The Supernova Cosmology Project" (SCP) composed 
by 580 type la supernovae [12j. The luminosity distance 
dh in a flat FRW cosmology, is computed through 



^(2, (0,(1, (2, #0) = c(l + z)H 



dz' 



£(z',Co,Ci,c 2 y 



FIG. 2: Confidence intervals for (("0,(2) of a bulk viscosity 
parametrization of the form C = Co + (2(0,/ a) (we set Ci — 0). 
We find as best estimates (Co = l-59±g;g|, C2 = 0.05lg; o f), 
see table U The solid red, dashed and filled contour plots 
correspond to the use of SNe, H(z) and the joint SNe + H(z) 
data sets respectively. The gray shaded area at the top left 
indicates the forbidden region where the total bulk viscosity 
Ctotai given by eq.© is negative at the present-day (z — 0). 
The confidence intervals shown correspond to 68.3%, 95.4% 
and 99.73% of probability. 



^(^,(0,(1, (2,^0) = m- M 



5 log 



in 



dL(zk, Co, (1, (2, Hq) 



Mpc 



+ 25, (15) 



where m and M are the apparent and absolute magni- 
tudes of the SNe la respectively, and the superscript 't' 
stands for "theoretical" . We construct the statistical \ 2 
function as 



XsNe(Ctb Cl) C2) = ^ 



[ fx* (z k , (0,(1,62) - Mfc 



(16) 



fc=i 



where \ik is the observational distance moduli for the fc-th 
supernova, a 2 is the variance of the measurement and n is 
the amount of supernova in the data set (n = 580). The 
results are shown in table U and the confidence intervals 
for the pairs (0 vs (1,(0 vs (2 and (4 vs (2 are shown in 
figures [T] to [3] respectively. The Hubble constant Hq is 
marginalized assuming a constant prior distribution (see 
appendix A of Q). 



where ^(0,(0,(17(2) is given by the expression ([8]) and 
c is the speed of light given in units of km/sec. The 
theoretical distance moduli for the fc-th supernova with 
redshift Zk is defined as 



B. Hubble expansion rate 

For the Hubble parameter H(z) measured at different 
redshifts, we use the 12 data listed in table 2 of Busca 



5 



Confidence Intervals 



1.5 



1.0 



0.5 



0.0 



-0.5 



■1.0 -, 




Case (<?(,,?,) 







FIG. 3: Confidence intervals for ((1,(2) of a bulk viscosity 
parametrization of the form ( = (l-H'+^a/a) (we set ( = 0). 
The solid red, dashed and filled contour plots correspond to 
the use of SNe, H(z) and the joint SNe + H{z) data sets 
respectively. The confidence intervals shown correspond to 
68.3%, 95.4% and 99.73% of probability. The best estimated 
values for ((1,(2) are all those points that lie on the line 
equation (2 = m(r+&, where m = -0.19808, b = 0.30918 with 
Xmin — 639.131 (green line) when using the combined SNe + 
H(z) data sets, except at the singular point ((1 = 1, (2 = 1/9) 
that corresponds to the vertex visible in the figure. The gray 
shaded area indicates the forbidden region where the total 
bulk viscosity (total given by eq.(© is negative at the present 
day(z = 0). 

et al. (2012) Q3, where 11 data come from references 
[ll-|l7|. The value H = 70 km s" 1 Mpc" 1 , is assumed 
for the data of Blake et al. (2011) [15| as Busca et al. 
suggest. The \ 2 function is defined as 
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xh(Co,Ci,C2) = 



g(z t ,Co,Ci,C 2 )-g° bs 

0~Hi 



(17) 



where H (z i: £ , Ci, C2) and H° hs are the theoretical and 
observed values respectively and ani the standard devi- 
ation of each H° hs entry. 

The total Xt function which combines the SNe and 
H{z) data sets together, is chosen in the following way: 



Xt 



XsNe 



where Xsnc ano - Xh are gi yen by expressions ([To]) and (fTT]) 
respectively. The function Xt i s then numerically min- 
imised in order to compute the "best estimates" for pairs 
of the viscous coefficients: (Co, Ci): ((0,(2), and (Ci,C 2 ), 
where the remaining viscous coefficient in each case is 
assumed to vanish. The minimum value of the x 2 func- 
tion gives the best estimated values of the pairs (Co,Ci)> 
(Co)C2)j and ((1,(2), and measures the goodness-of-fit of 
the model to data. 
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Case (f ,£j) 
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FIG. 4: Evolution of the total bulk viscosity ("total with re- 
spect to the redshift z, for the different cases. The plots 
come from the evaluation of the equation ([9]) at the best es- 
timated values for ((o,(i>C2) and error values. In the up- 
per and middle panels the red solid lines correspond to the 
evaluation of eq. © at ((0,(1,6) = (3.12,-0.96,0) and 
(1.59,0,0.05), respectively. The short dashed blue lines cor- 
respond to evaluation at (3.47,-0.66,0) and (1.65,0,0.057) 
respectively, while the long dashed black lines correspond to 
(2.78, -1.26, 0) and (1.54, 0, 0.04) respectively. In the bottom 
panel the best estimated values for ((1,(2), with (0 = 0, cor- 
respond to all those points which lie on the line (2 = m(i + b, 
where m = -0.19808, b = 0.30918 (see figure©. It turns 
out that the evolution of (total (z) is independent of the values 
of ((1,(2) when these are located at the best estimated line 
indicated above (the green line at figure©. The solid red line 
in the bottom panel corresponds to the evaluation of eq. @ 
at some arbitrary values of ((1,(2) related by the line equa- 
tion given above. The short dashed blue, and long dashed 
black lines correspond to the evaluation at the limits of the 
la confidence contour shown on figure© 



6 



Viscous model 



Best estimates 


Assumption 


2 2 
Xmin Xd.o.f. 


Co = 3.12±g-|| Ci = -0.96 ± 0.3 
Co = l-59±g : gg C2=0.05t° °o? 7 


C2 = 

Ci=o 


573.34 0.97 
573.34 0.97 



Case (£ ,£,) 



TABLE I: Best estimated values of the dimensionless viscous 
coefficients (Co,Ci,£2) for the dark component with bulk vis- 
cosity parametrized as C = Co + Ci-^ + ^ifl/a). It was as- 
sumed w = and Q v (z = 0) = 0.96 for the dark component, 
and flmo = 0.04 for the baryon matter. The best estimates 
were computed using the SNe + H(z) data sets. The first 
two columns show the best estimated values for pairs of vis- 
cous coefficients (Ci , Q) , where the remaining viscous coef- 
ficient Cfe in each case, is set to zero (third column). The 
fourth and fifth columns show the minimum of the \ 2 func- 
tion and its corresponding "x 2 function by degrees of free- 
dom": Xd.o.f. = Xmin/( n ~p)j where n is the number of data, 
and p is the number of free parameters estimated. The errors 
in the estimations are given to la. The best estimated values 
for (Ci, C2), with Co = 0, correspond to all those points which 
lie on the line C2 = m-Ci + b (m = -0.19808, b = 0.30918), with 
Xmin = 639.131 (xlo.f. = 1.08). Figures [1] to [3] show the con- 
fidence intervals. The Hubble constant Ho was marginalized 
assuming a constant prior distribution. 



The definition of "x 2 function by degrees of freedom''' 



Xd.o.f 



x 2 - , 

Amin/ 



/(n — p), where n is the number of total 
combined data used, and p is the number of free param- 
eters estimated, is also used in our computations. 




C. Local Second Law of Thermodynamics 

The local entropy production for a fluid on a FRW 
spacetime is expressed as [§] 

TV„s v = ({V„u u ) 2 = 9H 2 (, (18) 

where T is the temperature of the fluid, V ' v s v is the rate 
of entropy production in a unit volume, and £ is the total 
bulk viscosity. The second law of the thermodynamics 
can be stated as TV ' v s v > 0. Hence, since the Hub- 
ble parameter H is positive for an expanding Universe, 
Eq. (TT5|) implies that £ > 0, where £ is given by the ex- 
pression This inequality is an additional constraint 
in the possible values for the total dimensionless viscous 
parameter in our model. 



FIG. 5: Evolution of the dimensionless density parameters of 
the viscous dark matter component Q v (z) [solid red line] and 
of the baryon matter component Q, m (z) [dashed blue line] with 
respect to the redshift z, for the different cases. In the top and 
middle panels the curves correspond to the evaluation of eqs. 
(fT3|) at the best estimated values (Co, Ci, C2) = (3-12, —0.96, 0) 
and (1.59, 0, 0.05), respectively. In the bottom panel the best 
estimated values for (£i,Ca)j correspond to all those points 
which lie on the line: C2 = m£i +b, where m = —0.19808, b = 
0.30918 (see figure[3]). It turns out that the evolution of £l v (z) 
is also independent of the values of (Ci, C2) when these values 
are located at the best estimated line indicated above. The 
curves in the bottom panel correspond to the evaluation of 
eqs. (|13[) at some arbitrary values of (Ci, C2) related by the line 
equation given above. It is seen that starting at some redshift 
z ~ 7.7, fl„(z) > 1, while Sl v (z) < 0, which is a nonsense. 
From the top and middle panels it is seen that Q v (z) is always 
the dominant component with respect to Q m (z) at late times. 



IV. DYNAMICAL SYSTEMS 



The dynamical systems tools offer a powerful mean to 
extract relevant information out of the given cosmolog- 
ical model by investigating the equivalent phase space. 
Critical points in the phase space: past/future attrac- 
tors, saddle points, etc., can be correlated with generic 
solutions of the cosmological field equations. In order to 



be able to apply these tools one has to follow the steps 
enumerated here: i) to identify the phase space variables 
that allow writing the system of cosmological equations 
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FIG. 6: Evolution of the deceleration parameter q(z) with 
respect to the redshift z for the different cases evaluated 
at the best estimated values. In the top and middle pan- 
els the solid red lines correspond to the evaluation of eq. 
CE} at (Co, 6,6) = (3.12,-0.96,0) and (1.59,0,0.05), re- 
spectively, while the short dashed blue lines correspond to 
evaluation at (3.47,-0.66,0) and (1.65,0,0.057) respectively. 
The long dashed black lines correspond to (2.78, —1.26, 0) and 
(1.54,0,0.04) respectively. In the bottom panel the best es- 
timated values for (£i,C2), with £o = 0, correspond to all 
those points which lie on the line: C,2 = rnC,i + b, where 
m = -0.19808, b = 0.30918. The evolution of q(z) is in- 
dependent of the best estimated values of (Ci,C2). The solid 
red line in the bottom panel corresponds to the evaluation 
of eq. (|14p at some arbitrary values of (£i,m£i + b), while 
the short dashed blue and long dashed black lines correspond 
to the evaluation at the limits of the la confidence contour 
shown on figure [3] 



in the form of an autonomous system of ODE, 1 ii) with 
the help of the chosen phase space variables to build an 
autonomous system of ODE out of the original system of 
cosmological equations, and an usually forgotten or un- 
appreciated step, iii) to identify the phase space spanned 
by the chosen variables that is relevant to the cosmolog- 
ical model under study. 

Our goal in this section is to write the cosmological 
equations of the model <JSJ|4j) : 2 

Pm + 3Hp m = 0, 

p v + ZH{l + Lo) Pv -%H 2 = 0, 

3H 2 = p v + p m , 

6H + 6H 2 = -p m - (1 + 3uj) Pv + 9(H, (19) 

in the form of an autonomous system of ODE. To this 
end we have to choose appropriate phase space variables. 
In the present case our starting phase space variable is 
the dimensionless energy density parameter of the viscous 
fluid: 

x = tl v = < x < 1. (20) 

In terms of this variable the Friedmann constraint 
(third equation in (fT9"|) ) can be written as Q m = 1 — x, 
where we use the standard definition of the dimension- 
less energy density parameter of the i-th matter compo- 
nent, tti = pi/3H 2 . Also, one can write the following 
autonomous ODE: 

H' C 

x' = + lu)x - 2x— (21) 

or, since 

2^ = -3(l + u;z)+3|:, (22) 

the former equation can be written in more compact 
form: 



x' = 3(.t - 1) [ljx - jjj . (23) 

In the above equations the tilde accounts for derivative 
with respect to the parameter, r = In a. 

As already mentioned, here we shall investigate a vis- 
cous coefficient of the form given in equation ([S]): 

C = Co + Ciff + C 2 -, 

a 



1 There might be several different possible choices, however, not 
all of them allow for the minimum possible dimensionality of the 
phase space. 

2 Here we use units where 8ttG = 1. 
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or, since a/ a = H + H 2 , then the viscous parameter will 
obey the following equation: 



+ Cl + C2 + ~H C2 
= 2( /H + 2Ci - (1 + 3u;x)C 2 
2 - 3C 2 



(24) 



where, in the last row, we have taken into consideration 
equation ([22]). If we substitute back Eq.([24j) into (l23l) . 
we obtain the following master ODE: 



6 (a 



1) 



2-3C 2 



Co , . C2 



(25) 



Several cosmological parameters, such as the decelera- 
tion parameter q = — 1 — H' / H 1 and the equation of state 
(EoS) effective parameter u> e ff = — 1 — 2H' /3H, can also 
be rewritten in terms of the variable x. In fact, if take 
into account equations (|2"2")l and (IM1) one obtains: 



q = 



J eff 



I + 'Sojx- 3Co /H - 3(i 
2^3^ ' 
2ujx - 2Co/ff - 2Ci + C2 
2 - 3C 2 



(26) 



respectively. 

Depending on the case under consideration one would 
need yet another phase space variable which would be re- 
lated with the viscous coefficient Co (see below). In what 
follows we shall split the dynamical systems study into 
two different cases: A) when the viscous EoS parameter 
oj 7^ 0, and B) when the viscous fluid is dust: lo = 0. 



A. Viscous EoS w 7^ 

1. Case with C,o — 0. 

If we set Co = in Eg. (1231) . we obtain the following 
autonomous ODE for this particular case: 



, 6lo{x -I) ( 2Ci - C2 
x = — — — x - 



2-3C2 



2w 



(27) 



The phase space is the segment, ^ = {x\Q < x < \} . 
Two equilibrium/critical points are found. 



1. The first one, 

P v ■■ (a 



Pv 

3H 2 



1 , 



corresponds to the viscous matter-dominated solu- 
tion. The deceleration and EoS effective parame- 
ters, in this case, are given by: 

l + 3(w-Ci) 2(w-Ci) + C 2 

1 = n — 1 Ueff = 



The solution corresponds to accelerated expansion 
whenever, either £1 > (l+3w)/3, C2 < 2/3, or, £i < 
(1 + 3w)/3, C2 > 2/3. Otherwise it will correspond 
to decelerated expansion instead. 

If consider small r-dependent perturbation e = e(r) 
around this critical point: x — > 1 + e(t), up to 
C(e 2 ), the perturbation will obey the following lin- 
earized ODE: 

2w-2Ci + C2 , , 
e ( T > = 3 57 e ( r )' 

which can be readily integrated, 



e(r) = e e Ar , A = 3 



,C 2 -2Ci + 2^ 
2 - 3C 2 ' 



where eo is an integration constant. The solution 
is stable or, in other words, it is a future attractor 
in the phase segment if A < 0, i. e., if either 



2& - C2 > 2u, C2 < 2/3, 



2Ci - C2 < 2w, C2 > 2/3, 

which coincide with the regions in the space of pa- 
rameters (Ci, C2, where the point P v is correlated 
with inflationary expansion (q < 0). Otherwise, if 

2Ci - C2 < 2w, C2 < 2/3, 



2Ci - C2 > 2w, C2 > 2/3, 

the viscous matter-dominated solution is unstable 
(it is a past attractor), while the expansion occurs 
at a decelerated pace. Hence, either the viscous 
matter-dominated critical point P v (3H 2 — p v ), is 
the future inflationary attractor/end-point of any 
phase space orbit, or, alternatively, it is the past 
attractor/source point in the phase space, which is 
associated with decelerated expansion. 



2. The second critical point, 



m/v 



2Ci - C2 
2lo 



3H 2 



2u>p v 

2CT^C2' 



exists whenever 3 < 2£i — £2 < 2u, and corre- 
sponds to matter/viscous matter-scaling solution: 



a, 
a, 



2uj - 2Ci + C2 
2Ci - C2 



2-3C2 



2-3C2 



3 We shall be assuming that w is a non-negative quantity, which 
covers the most interesting physical situations. 
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The deceleration and the EoS parameters for this 
critical point are: q = 1/2 and w e // = respec- 
tively. 

Small r-dependent perturbations around this criti- 
cal point will obey, e' = Ae, or after integration, 



e(r) = e e x \ A = 



2Ci - C2 - 2u 
2 - 3C 2 



Hence the scaling critical point is stable, i. e., it is 
the future attractor in the phase space, if (2 < 2/3. 
Otherwise, if (2 > 2/3, this equilibrium point is 
the past attractor. In consequence, when the scal- 
ing equilibrium point exists, either it is the future 
attractor while the first equilibrium point (x = 1) 
is the past attractor, or vice versa. 

In this case (Co = 0), if 2£i — £ 2 < 2uj, then the orbits in 
the phase segment either depart from the viscous matter- 
dominated, decelerated solution, and end-up at the, also 
decelerated, scaling solution, or they are repelled from 
the scaling equilibrium point (corresponding to deceler- 
ated expansion always) and are attracted towards the 
viscous matter-dominated inflationary solution. In cither 
case none of these scenarios is suitable to accommodate 
the present cosmological paradigm, since there is no any 
critical point that could be associated with conventional 
matter and/or radiation dominance, which are included 
here in the matter component characterized by energy 
density /O m . 4 Radiation and matter-dominated phases 
are necessary to explain the formation of the amount of 
cosmic structure we see, in particular the right growth of 
structure (l3j . 

An interesting situation occurs when (2 — 0, Ci 0- 
In this case, for the viscous matter-dominated solution 
x = 1, one has 



1 + 3(w - Ci) 



J eff 



Ci 



so that the solution is inflationary if (1 > (1 + 3w)/3. Be- 
sides, this critical point is a future attractor in the phase 
space (segment) only if £1 > uj. Hence, this solution is 
stable and inflationary only if £1 > ui + 1/3. 

Alternatively, the scaling equilibrium point x = Qi/oj, 



^ - Ci 
Ci 



q = 1/2, u e ff = 0, 



exists whenever, < Ci < w - ft is stable whenever it 
exists, £1 < uj. Hence, when both critical points coexist, 
the viscous matter-dominated (decelerated) solution is 



4 As a matter of fact, in the present work, for simplicity, it has 
been assumed that the conventional matter behaves like pres- 
sureless dust, but it is clear that even if consider it to be ra- 
diation there would not be any critical point associated with 
radiation-domination (see the next section). 




FIG. 7: Existence of the critical points in parameters space 
(CiiCz) for the case IIV A 21 Regions in which the given(s) 
critical point (s) exists are the blue-coloured ones. The pa- 
rameter u has been arbitrarily chosen to be uj — 1/3 (viscous 
"radiation"). 



the past attractor, while the (also decelerating) scaling 
solution is the future attractor. Curiously, if Ci > uj + 
1, i. e., if the viscous matter-dominated solution is the 
future attractor (besides, it is the only critical point in 
the phase segment), the effective EoS parameter behaves 
like a phantom, uo e ff < — 1. 

If, on the contrary, (1 = 0, (2 7^ 0, then for the viscous 
matter-dominated critical point x — 1, 



q = 



1 + 3oj 2u + C2 



2-3C2' 



2-3C 2 



so that this solution corresponds to inflationary expan- 
sion if C2 > 2/3. It is stable when, either, (2 + 2uj < 0, 
C 2 < 2/3, or, (2 + 2w > 0, ( 2 > 2/3. The scaling solution, 
x = -C2/2W, 



2uj + ( 2 
C2 



exists if — 2lo < (2 < 0. It is stable whenever it exists. 
In this case the viscous matter-dominated cosmic fluid 
mimics phantom behavior if £2 > ^ + 1 (^e// < — 1). 

As before, no critical point in the phase space can 
be associated with conventional matter dominance as re- 
quired by the standard cosmological paradigm. 
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FIG. 8: In the upper panels the stability of the viscous matter- 
dominated critical point P v for case IIV A 21 is shown in the 
space of parameters (£i, £2) (as before £2 stands for the verti- 
cal axis while (1 is the horizontal axis, and uj = 1/3). The crit- 
ical point P v is stable in the blue-coloured regions in the upper 
left-hand figure, while, it is unstable in the blue-coloured re- 
gions in the upper right-hand figure. In the lower panel the 
regions in parameter space where the critical point P v cor- 
responds to inflationary phase are shown (also coloured in 
blue). 



2. Case with ( 0. 

In this case to the already existing variable x one has 
to add a new one: 



y 



1 Co = 1 - v 

Co/H+l H y 



Hence, the phase space is the bounded plane region 

* = {(x,y)\0 < x < 1, < y < 1}. 

As before, the Friedmann constraint reads Q TO = 1 — x. 
The corresponding autonomous system of ODE looks like 



*'= 3( !r 2 f*t C2 Vi) ( y + 2u}xy - 2 



(2 - 3C 2 ) y 
, _ 3 (2 - Ci - C2) 
2-3C 2 



y 



2 - 2Ci + C2 



For the deceleration and EoS parameters the following 
expressions are obtained: 



q 

Ueff 



(4-3Ci)j/ + 3^-3 
(2 - K2)y 
(2-2{ 1 +( 2 )y + 2cjxy-2 



(2 - 3C 2 )y 



(29) 



The critical points of (|28l) . Pi : (xi,yi), together with 
their main properties are summarized as follows. 

1. Viscous matter-dominated solution P v : (1,1). 
This case corresponds either to the formal limit 
Co = 0, or to the initial singular state character- 
ized by H — > 00. We have 



Q = 



1 - 3Ci + 3w 

2 - 3C 2 : 



J eff 



C2 - 2Ci + 2m 
2 - 3C 2 



The eigenvalues of the linearization (Jacobian) ma- 
trix for this point are: 



Ai = 
A 2 = 



3(1 -Ci-b+u) 

2 - 3C 2 
3(C 2 - 2Ci + 2uj) 
2 - 3C 2 



This solution is a past attractor if, either Ci < 
u> + 1/3, C 2 < 2/3, or Ci > w + 1/3, ( 2 > 2/3. 
The stability properties of this point in the space 
of parameters ( Ci : C2 ) , are shown in the figure 13 
where the free parameter w has been arbitrarily set 
equal to 1/3 ("viscous radiation"). 

2. de Sitter (also viscous fluid-dominated) solution 



P, 



ilS 



1 



H = Hn = 



2 - Ci - C 2 + w 
Co 



1 - Ci - C 2 + w 



characterized by q = —1, w e // = —1- The eigen- 
values of the Jacobian matrix corresponding to this 
critical point are: 



Ai = - 



3(1 -Ci -(2+L0) 

2 - 3C 2 
3 

(2-Ci-C 2 + w) 2 ' 



The de Sitter solution P^g exists whenever < y < 
1, i. e., if: Ci + C2 < ^ + 1 (see FigE]). It is the 
future attractor in ^ if: (1 + Q2 < u + 1, C 2 < 2/3. 
Otherwise, if: Ci + Ca < ^ + 1, £ 2 > 2/3, P dS 
is a saddle critical point instead. The stability of 
this equilibrium point in the space of parameters 
(Cij C 2 ) is shown in FigJHl where the additional free 
parameter lo has been arbitrarily chosen: lj = 1/3. 

3. Matter/viscous matter-scaling solution, 



p 



./V 



2C 



-C2 A On = 2uj - 2Ci + C 2 

2lo ' J ' n v 2Ci - C2 
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FIG. 9: Stability of the de Sitter critical point in the space 
of parameters ((1X2), for case IIV A 21 is shown in the upper 
panels while, in the lower panel the stability of the scaling 
point is depicted. The regions of interest in each case are the 
ones blue-coloured. 



This critical point exists (0 < x < 1) if (see FigEJ) , 

C2 , > , C2 

y <Ci<y+- 

It is characterized by, q = 1/2, v e ff — 0. The 
eigenvalues of the corresponding Jacobian matrix 
are 



Ai = 3/2, A 2 



3(( 2 - 2Ci + 2uj) 
2 - 3C 2 



so that, if, 2£i — C2 < 2w, C2 < 2/3, it is a saddle 
point in the phase space. Otherwise, if, 2£i — C2 < 
2cj, C2 > 2/3, it is the past attractor instead. In 
Fig El the stability of the scaling point is depicted in 
the space of parameters (Ci,C2)- As before viscous 
radiation has been considered (w = 1/3). 



B. Pressureless Viscous Fluid (oj — 0). 

In this case equations ([28)) simplify to: 

1 



a/= 3(2-26 +6) (lZ± ni , 



2-3C 2 
, 3(2 - d - C2) 



y 



Besides, 



2-3C 2 



y 



(y-i)U 



1-C1 + C2/2; ' 
1 



(4 - 3Ci)y - 3 
q= (2-3 C2 ) y >">" 



2 - Ci - C2 

(1 - Ci + (2/% - 1 
2(2 - 3C 2 )y 



(30) 



FIG. 10: Phase portraits (x,y) for the case IIV A 21 for differ- 
ent choices of the free parameters £1, and £2 (w = 1/3). The 
following values of the free parameters have been chosen in 
the upper panels: ft = — 1, £2 = — 2 - left-hand panel, and 
£1 = 2, C2 = 4 - right-hand panel. It is seen that, while in 
the latter case the de Sitter point does not exist (the past 
attractor/source point is the scaling solution and the viscous 
matter-dominated phase is a saddle point), in the former case, 
the viscous matter point is the past attractor, the scaling so- 
lution is a saddle critical point and the de Sitter phase is 
the future attractor. In the lower panels we have taken best 
estimated values of £1 = £i/3, and ("2 = 6^2 from tabled] 
Ci = -0.96, 6 = 0- left-hand panel, and Cl = 0, < 2 = 0.05 - 
right-hand panel. In both cases the orbits are repelled from 
the viscous matter-dominated solution and approach to the 
de Sitter future attractor. The scaling critical point does not 
exist. 



Only two critical points of the autonomous system of 
ODE ffl are found: 



1. Viscous matter-dominated solution P v : (1,1). The 
relevant parameters are: 

_ 1 - 3Ci ^ _ C2 - 2(i 
q 2 - 3C 2 ' Ueff 2 - 3C 2 ' 

while the eigenvalues of the linearization matrix 



3(C 2 -2Ci) 3(1 -Ci- 6) 

^1 — — ^ — 7 A 2 — 



2-3C 2 



2-3C 2 



The stability of this point in the space of parame- 
ters (Ci, C2), is shown in the figure ITU 



2. de Sitter equilibrium point 
1 



dS 



1, 



2 - Ci - Ca 



q = -1, UJ eff = -1, 
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FIG. 11: Stability of the viscous matter-dominated critical 
point in (£i, ^2)-space for the case IIV Bl (upper panels). In 
the lower panel the regions in the parameters space where 
this critical point is associated with inflationary expansion 
are shown. As before the regions of interest in each case are 
the blue-coloured ones. 



Fxiilenoe of de Sitter Point deSitten Poult Stable 




FIG. 12: Existence and stability of the de Sitter critical point 
in parameters space (£i, £2) for the case IIV Bl The regions in 
blue colour are the ones of interest in each case. 

which is also dominated by the bulk viscous matter 
(x = 1). The eigenvalues of the Jacobian matrix 
are: 

3(1 - Ci - C2) , 3 

2-3C 2 ' 2 (2-Ci -C2) 2 " 
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FIG. 13: Phase portraits (x, y) for the case IIV Bl for different 
choices of the free parameters (1, and £2 (w = 0). We have 
taken the best estimated values of £1 = Ci/3, and £2 = 6(^2 
(Co / 0) from table Q] Ci = -0.96, C2 = - left-hand panel, 
and fi = 0, (2 = 0.05 - right-hand panel. In both cases the 
orbits are repelled from the viscous matter-dominated solu- 
tion and approach to the de Sitter future attractor. 

The existence and stability of the de Sitter point in 
(C17 C2)-parameter space is depicted in Fig 021 

A distinguishing feature of the bulk viscous fluid sce- 
nario is the lack of equilibrium points in the phase 
space that could be correlated with conventional matter- 
dominance. A phase of conventional matter-dominance 
is required for the formation of the observed amount of 
cosmic structure. This is one of the most unwanted fea- 
tures of the viscous fluid scenario and, as shown, this 
conclusion is irrespective of whether we consider uj ^ 0, 
or, u> = 0. In the next section we shall see that this con- 
clusion is robust enough and it holds true even if add a 
radiation component to the model. 



V. MODEL OF BULK VISCOUS MATTER 
WITH CONVENTIONAL MATTER AND 
RADIATION 

In this section we shall investigate a more physically in- 
volved scenario with bulk viscous matter, with the bulk 
viscosity coefficient given by 1J5J. Here, besides a pres- 
sureless (ui — 0) viscous matter component, we shall in- 
clude conventional (non-relativistic or dust) matter, and 
also radiation. The cosmological equations are the fol- 
lowing: 

3H 2 = p r + p m + p v , 

6H + 6H 2 = -2p T - p m -p v + 9(H, 

p r + 4Hp r = 0, p m + 3Hp m = 0, 

Pv + 3H Pv - 9C# 2 = 0, (31) 

where p r is the energy density of the radiation component 
and, as before, p m and p v stand for the energy densities of 
non-relativistic (pressureless) matter and of bulk viscous 
(also pressureless) component, respectively. 
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In order to transform the above system of equations 
into a system of autonomous ODE we introduce the fol- 
lowing variables of the phase space (the two first variables 
x, and y, coincide with the former definitions and we add 
a new variable z): 



x = n v ,y 



1 



Pr 



(32) 



Co/h + v m 2 

The following autonomous system of ODE is obtained, 
x 1 = 3(1 — x )j^ + xz j 



z [ z- 1 



where 



H 



C_ = 2 + (2Ci-(z + l)C 2 -2)y 

H (2 - 3C 2 )y 



(33) 



(34) 



The Friedmann constraint can be written as, tt m = 1 — 
x — z, while the bounded 3D phase space is given by 

* = {(a;, y,z)\0<x<l, < y < 1, < z < 1}. (35) 

For the deceleration parameter q = —l — H'/H, and 
the effective EoS parameter 0J e ff — — 1 — 2H'/3H, one 
obtains, 



l + z-3(/H 



£_1 
3 #' 



respectively. 



A. Critical Points and Their Properties 

Three equilibrium points are found in the phase space 
^ (|35|) . A summary of these points, P$ : (xi,yi,Zi), to- 
gether with their main features is given below. 

1. Bulk viscous matter/radiation-scaling 

P r /v ■ (3(Ca - CO, 1, 1 - 3(Ca - CO) 

n r _ i - 3(c 2 - Ci) 
n, 3(c 2 - CO 

This solutions exists if 

o < 3(C 2 -CO < i, 

i. e., when Ci < C2 < Ci + 1/3- In this case 
the cosmic expansion is decelerating, q = 1, while 



u> e ff — 1/3. The eigenvalues of the Jacobian ma- 
trix for this point are, 



Ai — 1, A2 — 2, A3 — 



2[i-3(c 2 -CQ] 

2 - 3C 2 



so that it is a unstable critical point in ^ if C2 < 
2/3, and a saddle point otherwise. 

2. Bulk viscous matter-dominance, P v : (1,1,0) 
Cl v = 1. This point is characterized by 

l-3Ci C2 - 2Ci 

1 = o — > ^eff 



2 - 3C 2 



2 - 3C 2 ' 



and by the following eigenvalues of its linearization 
matrix: 

_ 2[1 - 3(C 2 ^ CO] x 3(1-Ci-Ca) 

^1 ^ ^ ) A 2 — 



A, 



2-3C2 
3(C 2 -2CQ 
2-3C 2 



2-3C 2 



Stability properties of this equilibrium point in the 
space of parameters (Ci, C2), are shown in figHH 

3. de Sitter (also bulk viscous matter-dominated, 
Q v — 1) solution, 



PdS : 1, 



1 



2 - Ci " C2 ■ 



, => H = 



Co 



1 - Ci - Ca 



For this solution q — <jJ e ff — — 1- Since the eigen- 
values of the Jacobian matrix are, 



Ai = -4, A 2 = -3, A 3 



3(1 ~ Ci ~ C2) 
2 - 3C 2 : 



then, whenever it exists, the de Sitter solution is 
a stable attractor in, Vf, if C2 < 2/3. Otherwise 
(C2 > 2/3), it is a saddle point instead (see FiglT4"|). 

A distinctive feature of this model is the absence of 
critical points associated with conventional dark matter 
and/or radiation domination. There is only one equi- 
librium point (P r / V ) where the radiation and the bulk 
viscous fluid scale in a constant fraction during the ex- 
pansion. 

As it is suggested by the results of the former sec- 
tion IIV1 if consider a different EoS parameter u ^ 
- say, tu = 1/3, - for the bulk viscous matter, then an 
additional conventional matter/bulk viscous fluid-scaling 
critical point P m / V would arise. Hence, in general, there 
can be found only critical points where the different com- 
ponents of the conventional matter (including radiation) 
scale with the viscous matter, but, in no case, conven- 
tional matter-dominance can be correlated with equilib- 
rium points in the phase space. This result can rule out 
the bulk viscous matter-dominated models as acceptable 
models for the description of the cosmological dynamics 
of our Universe. 
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FIG. 14: Parameters space (£1, £2) for the case studied in sec- 
tion[V] In the upper left-hand panel we show the region in the 
space of parameters where the bulk viscous matter/radiation- 
scaling solution exists (blue strip). Above the horizontal solid 
line at, Qi = 2/3, this solution represents a saddle point. 
Below this line it is a past attractor instead. In the up- 
per right-hand panel the blue region in (CijCa) is the one 
where the de Sitter solution exists. Above the horizontal 
solid line marking £2 = 2/3, the de Sitter solution is a sad- 
dle point. Below £2 = 2/3 it is a stable attractor. In the 
lower panels the regions in parameter space where the bulk 
viscous matter-dominated solution is either a stable attractor 
(left-hand panel), or a unstable point (right-hand panel), are 
coloured in blue. 

VI. DISCUSSION 

In this section we shall make use of the results of a 
combined investigation of type la supernovae test (sec- 
tion nO}, and of the dynamical systems study (sections 
IIVI andlV |) . to check the viability of the bulk viscous mat- 
ter model to explain the presently accepted cosmological 
paradigm. But, before going into the details of the dis- 
cussion, just a few comments on the role of the study of 
the equivalent phase space instead of the original set of 
cosmological equations. 

When we replace the original field variables, say H, 
p m , p v , etc., by the phase space variables x, y, etc., we 
have to keep in mind that, at the same time, we trade the 
original set of non-linear, second order differential equa- 
tions in respect to the cosmological time t (cosmological 
field equations), by a set of first order ordinary differen- 
tial equations with respect to the variable r = In a, 

x' = f(x,y,...), y' = g{x,y, ..),... 

The most important feature of the latter system of ODE 
is that the differential equations - alternatively, the func- 



tions f(x, y, ...), g{x,y ,...), etc. - do not depend explic- 
itly on the parameter r. We call the system autonomous. 
In other words, we are trading the study of the cosmolog- 
ical dynamics of H = H(t), p v = p v (t), etc., by the study 
of the flux in r-parameter of the equivalent autonomous 
system of ODE. Note that, if consider eternal expansion, 
< a < co, then t and r-chronologies are isomorphic; 

t S }to, 00 [ r G ] — 00, oo[ . 

The critical (also, equilibrium or fixed) points of this 
system, Pi : {xi,yi, ...), i. e., the roots of the system of 
algebraic equations 

f(x,y,...) = 0, g(x,y,...) =0,..., 

correspond to (classes of) solutions of the original sys- 
tem of cosmological equations. If a given equilibrium 
point P a : ( ) is a (global) stable attractor, i. e., 

small linear perturbations around P a , x — > x a + 
V Ua + ({t), etc., decay with r, 5 then, independent on 
the initial conditions chosen, x(tq) = xq, y(ro) — y$,..., 
every orbit in the phase space will approach P a into the 
future ({r : r > t }), i. e., the stable (global) attractor 
is the end point of any orbit in \E'. On the contrary, if a 
given critical point, P s : (x Sl y s , ...), were unstable, i. e., 
small perturbations around, P s , uncontrollably grow up 
with r, then this point were a past attractor or, also, the 
source point of any orbit in the phase space. For a third 
class of critical points; the so called "saddle points", de- 
pending on the initial conditions chosen, orbits in ^ can 
approach to this point, spend some time around it, and 
then be repelled from it to finally approach to the stable 
attractor (if it exists). 6 

Now suppose we have a typical phase portrait, com- 
posed of a source critical point P s , a saddle point P* , and 
a stable (global) attractor P a . Each one of these points 
corresponds to given solutions of the original cosmologi- 
cal equations, H — H s (z), H — H*(z), and H — H a (z), 
respectively. In the above expressions z is the redshift, 
which is related with the parameter r through 

1 - a(t) 

t = — ]n(z + 1), z = z{t) = — 



5 This occurs if the eigenvalues of the Jacobian matrix 

(h fv -\ 

V ■ '■■ ..J 

all have negative real parts. 

6 For purposes of space, our discussion here is oversimplified, since, 
in general, critical points can be of many types, for instance, 
spiral, etc. Besides, there can be found also (un)stable manifolds 
such as cycles, etc. To worsen things there can coexist several 
local attractors, saddle points, etc., so that, in general, a given 
orbit in the phase space can approach to several saddle points 
before they end up at a given local attractor. 
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A also typical orbit in the phase space will start at P s 
for t = — oo, then will approach to P*, and after a finite 
(perhaps sufficiently long) At, will be repelled by P* to 
finally be attracted towards P a . The parallel history in 
terms of the equivalent cosmological dynamics will be the 
following. The expansion starts with a Hubble parame- 
ter dynamics H = H s (z), then, as the Universe expands, 
the cosmic history enters a transient period character- 
ized by the dynamics dictated by H = P*(z). After a 
perhaps long yet finite period Az, the cosmic expansion 
will abandon the latter phase, to enter into a stage which 
dynamics obeys H = H a (z) lasting for ever. 

To illustrate the above discussion with a concrete 
cosmological model, take as an example the so called 
"lambda-cold-dark-matter" (A-CDM) model. The cos- 
mological equations for this model in a flat FRW space- 
time are, 

3P 2 = p m + A, 

2H = -p m , p m + 3Hp m = 0, 

where p m is the energy density of the CDM and A is 
the cosmological constant. It is convenient to introduce 
the variable x = A/3P 2 of the phase segment {x\0 < 
x < 1}. The Friedmann constraint can be written as 
fi TO = 1 — x. The autonomous ODE obtained in this case 
is the following: x' = 3x(l — x). There are two critical 
points of this ODE: i) the matter dominated solution 
x = (fi m = 1), which happens to be a source point 
(past attractor), and ii) the de Sitter solution x = 1 (H = 
-\/A/3), which is the stable (future) attractor. Hence, 
since a typical trajectory in the phase segment starts at 
x = and end ups at x = 1, then in this model the cosmic 
history starts in a matter dominated period (essential for 
the formation of structure) and end ups in an inflationary 
stage lasting for ever into the future. 

If we improve this model by adding a radiation mat- 
ter component, there would be three critical points in 
the (now 2D) phase space: i) a radiation-dominated 
phase (unstable critical point), followed by ii) a tran- 
sient matter-dominated stage (saddle point), and iii) a 
de Sitter point (stable attractor). In this case the cosmic 
history starts at a radiation-dominated stage, then enters 
a transient period of matter dominance (essential for the 
formation of cosmic structure) , to finally approach to the 
stable attractor (the de Sitter phase) which will last for 
ever into the future. This is, precisely, the behavior one 
expects from a model designed to recreate the presently 
adopted cosmological paradigm. One has to care only 
about giving appropriate initial conditions xq = x(tq), 
etc., so as to get enough formation of structure, i. e., to 
ensure that the corresponding orbit in the phase space 
will spend enough time in the neighbourhood of the sad- 
dle critical point associated with the matter-dominated 
solution. 

We want to underline that, the fact that a given exact 
solution of the cosmological equations can not be associ- 
ated with a critical point in the equivalent phase space, 



does not mean at all that the above solution does not 
exist. In fact it might exist, but it should not be a sta- 
ble solution, so it should not be useful in a cosmologi- 
cal setting. To illustrate our point take as an example, 
again, the A-CDM model. There is a large set of known 
(classes of) solutions to the cosmological equations of the 
model, however, only the radiation-dominated, matter- 
dominated, and the de Sitter solutions, are of importance 
in a cosmological context which is compatible with the 
accepted cosmological paradigm. The latter is charac- 
terized by the following stages: i) early time inflationary 
period, followed by ii) a radiation-dominated, and iii) a 
matter-dominated phases, both associated with a stage 
of decelerated expansion, and iv) a present period of ac- 
celerated expansion which might last forever. 

What happens in the model of interest in this paper? 
Here, without loss of generality, we shall analyse only 
the model studied in section [V] where the cosmological 
dynamics is fuelled by radiation, conventional pressure- 
less matter, and bulk viscous (also pressureless) mat- 
ter. There are found three critical points (Co ^ 0) of 
the equivalent autonomous system of ODE correspond- 
ing to this model: i) bulk viscous matter/radiation- 
scaling solution P v / r , ii) bulk viscous matter dominance 
P v , and iii) de Sitter (also bulk viscous matter domi- 
nated) solution Pds- If set Co = only the bulk vis- 
cous matter/radiation-scaling solution (point P v / r ), and 
the bulk viscous matter-dominated point P v , survive. In 
what follows we shall focus in the physically more inter- 
esting case where Co 7^ which shows a richer phase space 
dynamics. Regions in the space of parameters (Ci , C2) ~ 
case with Co 7^ _ where the above critical points ex- 
ist as well as their stability properties are shown in the 
figure [HI It is seen that while the point P v / r exists in 
a very narrow strip in the phase space, the bulk viscous 
matter-dominated point P v always exists. 

There are four regions of interest in the parameters 
space (Ci, C2): 

• Rl: region where the three points co-exist to- 
gether. In this region the viscous matter/radiation- 
scaling solution (P v / r ) is the past attractor, the 
bulk viscous matter-dominated phase (P v ) is a sad- 
dle point, while the de Sitter solution (Pds) is the 
stable attractor. 

• R2: region where only P v / r and P v co-exist. In 
this case the bulk viscous matter-dominated solu- 
tion (point P v ) is the past attractor, while, the vis- 
cous matter/radiation-scaling solution is a saddle 
point in the phase space. Not of cosmological in- 
terest since there is not critical point which can be 
associated with a present stage of accelerated ex- 
pansion. 

• R3: region where only the critical point P v - corre- 
sponding to viscous matter dominance - and the de 
Sitter solution (point Pds) co-exist. In this case, ei- 
ther i) P v is the past attractor and Pds is the stable 
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attractor, or ii) P v is a saddle point while Pds con- 
tinues being the future attractor (no past attractor 
exists), or iii) Pds is a saddle point while P v is the 
stable attractor (no past attractor exist). The lat- 
ter case is not of interest for cosmology since there 
is not any critical point which could be associated 
with decelerated expansion. 

• R4: region where the only existing critical point is 
P v - bulk viscous matter dominance. In this region 
P v is either a saddle or an unstable critical point. 
Not of interest for cosmology since there is not any 
critical point corresponding to the present stage of 
accelerated expansion. 

As seen, only regions Rl, R3 i), and R3 ii) could be 
of cosmological interest. If substitute £i and £2 by their 
best estimated values in table HI in both cases the region 
of parameters space Rl above is single out. In this case 
Py/ r is the past attractor, P v is a saddle point (both P v / r 
and P v are associated with decelerated expansion), while 
PdS is the stable attractor. As mentioned, in Rl there 
is an equilibrium point (point P v / r ) where the bulk vis- 
cous fluid scales in a constant fraction with radiation. In 
this case the effective fluid behaves like radiation, but one 
which is partly bulk viscous. As already mentioned in the 
introduction, in a cosmological setting the bulk viscosity 
may arise when the cosmic fluid expands (or contracts) 
too fast so that the system does not have enough time to 
restore its local thermodynamic equilibrium and, then, 
it arises an effective pressure restoring the system to its 
thermal equilibrium. When the fluid reaches again the 
thermal equilibrium the bulk viscous pressure vanishes 
0, @]. This would mean that the stage of the expan- 
sion associated with P v / r may last for just a brief period 
of time while the Universe was out of thermodynamic 
equilibrium, perhaps not enough to produce the right 
peace of the growth of the fluctuations. If there is a very 
brief period of radiation domination, then, as long as the 
fluctuations re-enter the horizon these will not be suf- 
ficiently damped, and correspondingly, an unacceptable 
large value of the dispersion of the density contrast at 
the scale 8 h~ 1 Mpc (cr$) might be obtained (see a similar 
discussion but in the opposite direction in Ref.[13]). Be- 
sides, in either case although in the bulk viscous matter- 
dominated stage the effective fluid mimics dust, viscosity 
may affect the formation of structure in a way that can 
be observationally tested. 

As seen from the above analysis the most serious ob- 
jection against bulk viscous matter scenarios is the ab- 
sence of conventional matter and radiation dominated 
eras. Such a behavior is in marked contradiction with 
the big bang paradigm according to which, back enough 
into the past when the temperature of the Universe was 
larger than 10 4 K, the dynamics of the cosmic evolu- 
tion was driven by a relativistic mixture in the form 
of radiation (radiation dominated stage). As the Uni- 
verse expanded and cooled down to temperatures below 



10 4 K (and up to 3 x 10 3 K), the density of radiation 
diluted and the cosmic evolution entered in a stage of 
(non-relativistic) matter dominance. During this phase 
radiation decoupled from baryons to form the cosmic mi- 
crowave background and, what is more important, atoms 
and the derived cosmic structure we see (galaxies, clus- 
ters of galaxies, etc.) were formed. 

Here we have shown that the absence of such conven- 
tional matter and radiation-dominated phases is irrespec- 
tive of the region in the space of parameters Co, Cii C2, 
and lu, so that our conclusion is robust enough and the 
bulk viscous matter-dominated scenario should be ruled 
out, at least for the parametrization considered in this 
paper. 



VII. CONCLUSION 

In this paper we have applied the dynamical systems 
tools, in conjunction with the SNe la data testing, to 
judge about the possibility that cosmological bulk viscous 
matter can stand for an alternative to dark energy. We 
chose a formerly used parametrization of bulk viscosity 
0, [1] with the addition of a term measuring the influence 
of the acceleration of expansion: 



C - Co + CiH + C2 




The study of the asymptotic properties of the model 
in the equivalent phase space shows that there are 
not critical points that could be associated with either 
conventional radiation or matter dominance. This re- 
sult is independent of the values taken by the free pa- 
rameters of the model. In consequence, the bulk vis- 
cous matter-dominated model is not able to accom- 
modate the presently accepted cosmological paradigm. 
This argument alone can be considered as a serious 
objection against cosmological models of bulk viscos- 
ity. Notwithstanding, we recommend that other possible 
parametrizations of bulk viscosity should be considered 
before concluding to rule out the models. 
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